Monopole giant resonances and nuclear 
compressibility in relativistic mean field 

theory 



D. Vretenar 

Physics Department, Faculty of Science, University of Zagreb, 

Zagreb, Croatia 

G. A. Lalazissis 

Department of Theoretical Physics, Aristotle University of Thessaloniki 

GR-54006 Thessaloniki, Greece 

R. Behnsch, W. Poschl and P. Ring 
Physik-Department der Technischen Universitat Miinchen, 
D-85748 Garching, Germany 

February 9, 2008 



Abstract 

Isoscalar and isovector monopole oscillations that correspond to giant resonances 
in spherical nuclei are described in the framework of time-dependent relativistic 
mean-field (RMF) theory. Excitation energies and the structure of eigenmodes are 
determined from a Fourier analysis of dynamical monopole moments and densities. 
The generator coordinate method, with generating functions that are solutions of 
constrained RMF calculations, is also used to calculate excitation energies and tran- 
sition densities of giant monopole states. Calculations are performed with effective 
interactions which differ in their prediction of the nuclear matter compression modu- 
lus K-ara- Both time-dependent and constrained RMF results indicate that empirical 
GMR energies are best reproduced by an effective force with !^ 270 MeV. 
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1 Introduction 



The nuclear matter compression modulus i^nm is an important quantity in the description 
of properties of nuclei, supernovae explosions, neutron stars, and heavy ion collisions. In 
principle the value of -ft'nm can be extracted from the experimental energies of isoscalar 
monopole vibrations (giant monopole resonances GMR) in nuclei. In a semi empirical 
macroscopic approach, the analysis is based upon a leptodermous expansion of the com- 
pression modulus of a nucleus, analogous to the liquid drop mass formula [|I| . In principle 
such an expansion provides a "model independent" determination of K^^. However, the 
macroscopic approach presents several ambiguities. The formula itself is based on the 
assumption that the breathing mode is a small amplitude vibration. More important, to 
correctly interpret the value of the volume term in the expansion, one has to assume a 
certain mode of vibration. It has been argued that a direct determination of the various 
contributions to the compression modulus through a fit of the breathing mode frequencies 
cannot provide an accurate value for -ft'nm- Recently it was shown by Shlomo and Young- 
blood that the complete experimental data set on isoscalar monopole giant resonances 
does not limit the range of -R'nm to better than 200 — 350 MeV. Microscopic calculations 
of GMR excitation energies might provide a more reliable approach to the determination 
of the nuclear matter compression modulus. One starts by constructing sets of effective 
interactions which differ mostly by their prediction for i^nm, but otherwise reproduce 
reasonably well experimental data on nuclear properties. Static properties of the ground 
state are calculated in the self-consistent mean field approximation, and RPA calculations 
are performed for the excitations P, ^. Non relativistic Hartree-Fock plus RPA calcula- 
tions, using both Skyrme and Gogny effective interactions, indicate that the value of K^ra 
should be in the range 210-220 MeV. The breathing-mode giant monopole resonances 
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have also been studied within the framework of the relativistic mean field (RMF) theory, 
using the generator coordinate method (GCM) [§. It was shown that the GCM succeeds 
in describing the GMR energies in nuclei and that the empirical breathing mode energies 
of heavy nuclei can be reproduced by effective forces with i^Tnm ~ 300 MeV in the RMF 
theory. 

In the present article we describe isoscalar and isovector monopole oscillations in 
spherical nuclei in the framework of time dependent relativistic mean field theory (RMFT). 
The model represents a relativistic generalization of the time-dependent Hartree-Fock ap- 
proach. Nuclear dynamics is described by the simultaneous evolution of A single particle 
wave-functions in the time-dependent mean fields. Frequencies of eigenmodes are deter- 
mined from a Fourier analysis of dynamical quantities. In this microscopic description, 
self-consistent mean-field calculations are performed for static ground-state properties, 
and time-dependent calculations for monopole excitations using the same parameter sets 
of the Lagrangian. A basic advantage of the time dependent model is that no assumption 
about the nature of the mode of vibrations has to be made. Another approach that goes 
beyond the HF-I-RPA approximation is provided by the Generator Coordinate Method. 
In the second part of the article we extend the model of Ref . , and calculate the exci- 
tation energies of giant monopole states with relativistic GCM. As generating functions 
we use products of Slater determinants , built from single-particle solutions of constrained 
RMF calculations, and coherent states that represent the meson fields. 

The article is organized as follows. In Sec. 2 we present the essential features of 
the time-dependent relativistic mean-field model, and some details of its application to 
spherical nuclei. Results of time-dependent calculations for a series of doubly closed-shell 
nuclei and a set of effective interactions are discussed in Sec. 3. In Sec. 4 details of the 
relativistic GCM are explained. Sec. 5 contains a discussion of the systematics of energies 
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of giant monopole states that result from constrained RMF calculations with GCM. A 
summary of our results is presented in Sec. 6. 

2 The time-dependent relativistic mean-field model 

The dynamics of collective motion in nuclei is described in the framework of relativis- 
tic mean- field theory ^ |, |I^. Details of the time- dependent model are given in 
Refs. [0, Here we only outline its essential features. In relativistic quantum hadro- 
dynamics the nucleons, described as Dirac particles, are coupled to exchange mesons and 
photons through an effective Lagrangian. The model is based on the one boson exchange 
description of the nucleon-nucleon interaction. The Lagrangian density of the model is 
given as 



- gaixJip - gco-ipl ■ uJip - gpip'j ■ prip - eip'j ■ A^^^—^^ip . (1) 

The Dirac spinor if) denotes the nucleon with mass m. m^, m^, and rup are the masses 
of the cr-meson, the cj-meson, and the p- meson, and g^, g^, and gp are the corresponding 
coupling constants for the mesons to the nucleon. U{a) denotes the nonlinear a self- 
interaction 

U{a) = \mla' + \g,a' + \g,a\ (2) 

and nt"", R''^, and F''" are field tensors 0. 

From the Lagrangian density (|1|) the coupled equations of motion are derived, the 
Dirac equation for the nucleons: 



idtipi 



ct l-iV - g^uj - gpTp - e^^^—^A] + P{m + g„a) 



(1 - T-s) / 
+gu^^o + QpTpQ + e Aq 



(3) 



and the Klein-Gordon equations for the mesons: 

(d^-A + ml)a = -g^ps - g2a^ - gscr^ (4) 
[d^ - A + ml)uj^ = gj^ (5) 



dt -A + rrip) = gpj^ (6) 
(5' -a) A, = ej^. (7) 

In the relativistic mean field approximation the A nucleons, described by a Slater determi- 
nant 1$) of single-particle spinors ipi, [i = 1,2, A), move independently in the classical 



meson fields. The sources of the fields are calculated in the no-sea approximation |12 
- the scalar density 



the isoscalar baryon current 



the isovector baryon current 



A 

ps = X^V'iV'i, (8) 

i=l 
A 

f = E^iY^i, (9) 

i=l 
A 

= E^^T'^r^i, (10) 

i=l 



- the electromagnetic current for the photon-field 

Ji^m = (H) 
i=l ^ 

where the summation is over all occupied states in the Slater determinant |$). Negative- 
energy states do not contribute to the densities in the no-sea approximation of the sta- 
tionary solutions. However, negative energy contributions are included implicitly in the 
time-dependent calculation, since the Dirac equation is solved at each step in time for a 
different basis set. Negative energy components with respect to the original ground state 



basis are taken into account automatically, even if at each time the no-sea approxima- 
tion is applied. It is also assumed that nucleon single-particle states do not mix isospin. 
Because of charge conservation, only the 3-component of the isovector p contributes. 

The ground state of a nucleus is described by the stationary solution of the cou- 
pled system of equations It specifies part of the initial conditions for the time- 
dependent problem. In the present work we consider doubly closed-shell nuclei, i.e. nuclei 
that are spherical in the ground state. The nucleon spinor is in this case characterized by 
the angular momentum j, its 2;-projection m, the parity vr, and the isospin ^3 = ±| for 
neutrons and protons [^] 

^h(v tsU) = -( ^(^)*'^-^^' ^' ^ e-'^' y (U) (12) 

Xt is the isospin function, the orbital angular momenta / and / are determined by j and 
the parity vr, /(r) and g{r) are radial functions, and ^ijm is the tensor product of the 
orbital and spin functions 

<S>ijm{9,^,s) = <^ms I mi\j m>Yim,{9,ip) Xms{s). (13) 

nismi 

For a given set of initial conditions, i.e. initial values for the densities and currents 
in Eqs. (P|-pT]), the model describes the time evolution of A single particle wave-functions 
in the time-dependent mean fields.. The description of nuclear dynamics as a time- 
dependent initial- value problem is intrinsically semi-classical, since there is no systematic 
procedure to derive the initial conditions that characterize the motion of a specific mode 
of the nuclear system. In principle the theory is quantized by the requirement that there 
exist time-periodic solutions of the equations of motion, which give integer multiples of 
Planck's constant for the classical action along one period [|13|. For giant resonances the 
time-dependence of collective dynamical quantities is actually not periodic, since generally 
giant resonances are not stationary states of the mean-field Hamiltonian. The coupling 
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of the mean-field to the particle continuum allows for the decay of giant resonances by 
direct escape of particles. In the limit of small amplitude oscillations, however, the en- 
ergy obtained from the frequency of the oscillation coincides with the excitation energy of 
the collective state. In Refs. |jl2|, we have studied oscillations of isovector dipole. 



isovector spin-dipole, isoscalar quadrupole, and isovector quadrupole character, which cor- 



respond to giant resonances in spherical nuclei. In Ref. |T3| double giant isovector dipole 
and double giant isoscalar quadrupole resonances have been described in the framework 
of the time-dependent RMF theory and compared with recent experimental results. The 
model reproduces reasonably well the experimental data on energies and, for light nuclei, 
the widths of giant resonances. All the results are obtained using parameter sets that 
reproduce ground-state properties of nuclei, i.e. no new parameters are introduced in the 
model to specifically describe giant resonances. In the present study we apply the model 
to isoscalar and isovector monopole oscillations in spherical nuclei. In this microscopic 
description, self-consistent mean-field calculations are performed for static ground-state 
properties, and time- dependent calculations for monopole excitations. Using sets of effec- 
tive interactions which differ mostly by their prediction of nuclear matter compressibility, 
and which otherwise provide reasonable results for ground state properties, we calculate 
the excitation energies of monopole resonances in a series of spherical nuclei. 

In order to excite monopole oscillations in a doubly closed-shell nucleus, the spherical 
solution for the ground-state has to be initially compressed or radially expanded by scaling 
the radial coordinate. The amplitudes and g'^"'^ of the Dirac spinor are defined 

r°"(r^on) = (1 + a) /(r), <7"^°'^(r^on) = (1 + a) g{r) . (14) 

The new coordinates are 

rmon = (1 + a) r. (15) 
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A reasonable choice for the parameter is \a\ ~ 0.05 — 0.1. The energy transferred in this 
way to the ground state of the nucleus is somewhat above the giant resonance energy, 
and the resulting oscillations do not show anharmonicities associated with large amplitude 
motion. For the case of isoscalar oscillations the monopole deformations of the proton and 
neutron densities have the same sign. To excite isovector oscillations, the initial monopole 
deformation parameters of protons and neutrons must have opposite signs. After the 
initial monopole deformation (|l^, the proton and neutron densities have to be normalized. 
It should also be emphasized that apart from the fact that we concentrate on monopole 
excitations no assumption about the radial nature of the mode of vibrations is made in 
the time-dependent calculation. We do not have to assume that the motion is adiabatic, 
nor that the mode corresponds to a scaling of the density. The frequency dependence of 
dynamical quantities and the transition densities are used to determine the structure of 
the eigenmodes. The frequency can be simply related to nuclear compressibility only if a 
single compression mode dominates. 

For the case of spherical symmetry, the time-dependent Dirac equation (j^) reduces 
to a set of coupled first-order partial differential equations for the complex amplitudes / 
and g of proton and neutron states 

^dtf = {Vo + g.cr)f+{dr---iVr)g (16) 

r 

idtg = {Vo - ga(J - '2m)g - {dr + ^ - iK)/, (17) 

where k = ±(j + |) for j = / =F |, and the indices and r denote the time and radial 
components of the vector field 

= g^u^ + g,p^,T, + e^^^A^ (18) 
For a given set of initial conditions, the equations of motion propagate the nuclear system 
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in time. The potentials are solutions of the Klein-Gordon equations 

( 2 d 



+ m^U(r) = s^(r), (19) 



are meson masses for (p = a, u and p, and zero for the photon. The source terms are 
calculated from ®-(|ll|) using in each time step the latest values for the nucleon ampli- 
tudes. Retardation effects for the meson fields are not included, i.e. the time derivatives 
in the equations of motions for the meson fields are neglected. This is justified by the 
large masses in the meson propagators causing a short range of the corresponding meson 
exchange forces. At each step in time the meson fields and electromagnetic potentials are 
calculated from 

POO 

0(r) = / G^ir,r')s^iry''dr' , (20) 
Jo 

where for massive fields 

G^r,r') = -^i^(e-^l-'-'l_e-.l'^+^'-'l), (21) 
2m^ rr ^ ' 

and for the Coulomb field 

Gc{r, r') = - for r > r' 
r 

Gc{r, r') = for r < r'. (22) 

The dynamical variables that characterize vibrations of a nucleus are defined as expec- 
tation values of single-particle operators in the time-dependent Slater determinant |$(t)) 
of occupied states. In the framework of the TDRMF model the wave-function of the 
nuclear system is a Slater determinant at all times. For isoscalar monopole vibrations, 
the time-dependent monopole moment is simply defined as: 

(At)) = i($(t)|r^|$(t)). (23) 



The corresponding Fourier transform determines the frequencies of eigenmodes. The tran- 
sition density for a specific mode of oscillations is defined as the Fourier transform of the 
time-dependent baryon density at the corresponding frequency. The effective compression 
modulus Ka for a nucleus of mass number A is defined as 



where Eq is the energy of the isoscalar giant monopole resonance, m is the nucleon mass, 
and < >o denotes the expectation value of in the unperturbed ground state. As 
in Ref. 0], Eq. (l2^ ) represents just a convenient parameterization of the dominant A- 
dependence for the energy of the giant monopole state. 

3 Monopole oscillations in spherical nuclei 

The study of isoscalar monopole resonances in nuclei provides an important source of in- 
formation on the nuclear matter compressibility. We first consider three spherical nuclei 
from different regions of the periodic table: ^°Ca, ^°Zr and ^°^Pb. These nuclei differ not 
only in their masses, but also in the ratio of proton to neutron number. The complete 
experimental data set on isoscalar monopole giant resonances has been recently analyzed 
by Shlomo and Youngblood p. Using a semi empirical macroscopic approach to deduce 
from the systematics of GMR the systematics of nuclear compressibility, and to extrapo- 
late the data to infinite nuclear matter, they have shown that the complete data set does 
not limit the range of K^^ to better than 200 - 350 MeV. The GMR energy in ^ospb is 
rather well established at 13.7 ± 0.3 MeV. The average GMR energy for ^°Zr as deduced 
from various experiments is around 16.5 MeV, and for ^''Ca the value of the excitation 
energy adopted in the calculation was 18 ± 1 MeV [0. In light nuclei, of course, not all 
monopole strength is seen. 



Ti^AKa 



(24) 
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We have performed time-dependent relativistic mean-field calculations for six sets 



of Lagrangian parameters: NLl NL3 NL-SH |3, NL2 IT8|, HS and LI [|g] 
(in order of increasing value of the nuclear matter compression modulus -ft'nm)- The values 
of the parameters are given in Table 1 . These effective forces have been frequently used to 
calculate properties of nuclear matter and of finite nuclei. In this respect our work parallels 
the generator coordinate calculations of giant monopole resonances and the constrained 
compressibility within the relativistic mean- field model of Ref . , and the non- relativistic 
self-consistent mean-field plus RPA calculation of nuclear compressibility with Gogny 
effective interactions of Ref. 0] . The idea is to restrict the possible values of the nuclear 
matter compression modulus, on the basis of the excitation energies of giant monopole 
states calculated with different effective interactions. In addition to the four non-linear 
sets NLl, NL3, NL-SH and NL2, we also include two older linear parameterizations HS 
and LI, which do not include the self-coupling of the sigma field. The sets NLl, NL-SH 
and NL2 have been extensively used in the description of properties of finite nuclei |^ [l^ . 
In order to bridge the gap between NLl {K^^ = 211.7 MeV), and NL-SH {K^^ = 355.0 
MeV), we have also included a new effective interaction NL3 (-ft'nm = 271.8 MeV). This 
new parameter set has been derived recently |]T6l by fitting ground state properties of a 
large number of spherical nuclei, as well as nuclear matter properties. It appears that the 
NL3 effective interaction reproduces experimental data better than, for instance, NLl or 
NL-SH. 

Results of model calculations for ^°Ca are displayed in Figs, la and lb: time- 
dependent monopole moments and the corresponding Fourier power spectra are shown. 
Time is measured in units of fm/c, and the energy E = hu in MeV. The numerical ac- 
curacy is AE = 27r/ic/Tfinai = 27r^c/1000fm ^ 1.2 MeV. The frequency of the isoscalar 
monopole oscillations generally increases with the nuclear matter compression modulus. 
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and the vibrations become more anharmonic. For the first three sets of parameters (NLl, 
NL3, and NL-SH) a single mode dominates. For NL2, HS and LI the Fourier spectra 
are fragmented. The NLl effective interaction reproduces especially well the ground-state 
properties of nuclei close to the stability line 0, and it appears that the frequency of 
monopole oscillations for the NLl parameter set is very close to the expected experimen- 
tal excitation energy of 18 ± 1 MeV. 

In Figs. 2a and 2b the results for ^°Zr are shown. The frequencies increase with 
i^'nm; the motlou is more anharmonic and more damped. The damping of the monopole 
moments comes from the coupling to the continuum, and there is also a contribution from 
the damping via the mean- field (Landau damping). The experimental GMR excitation 
energy for ^^Zt of 16.5 MeV is found between the values calculated for the NLl and 
NL3 parameter sets. For NLl (15.6 MeV), NL3 (19.4 MeV), NL-SH (21.1 MeV) and 
the lower frequency of NL2 (21.6 MeV) we display in Fig. 3 the corresponding transition 
densities. For NLl, NL3 and NL-SH the transition densities show a radial dependence 
characteristic for the "breathing" mode: a change of the density in the volume at the 
expense of that on the surface. The details depend on the underlying shell structure. An 
unusual radial dependence is found for the NL2 parameterization. The transition density, 
both for protons and neutrons, has a minimum at r = 0, of the same sign as on the 
surface. The monopole oscillations at the center of the nucleus are in phase with surface 
oscillations. It seems that a similar behavior is found for HS and LI, although due to 
fragmentation it is difficult to decide at which frequency to plot the transition density. 

In Ref. 0] it was stressed that, rather than on the systematics over the whole periodic 
table, the determination of the nuclear compressibility relies more on a good measurement 
and microscopic calculations of GMR in a single heavy nucleus such as ^^'^Pb. The results 
of TDRMF calculations for ^osp]^ presented in Figs. 4-7. In Figs. 4a and 4b we 
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display the monopole moments and their Fourier power spectra. As one would expect for 
a heavy nucleus, there is very little fragmentation and a single mode dominates for all 
parameter sets. The experimental excitation energy 13.7 ± 0.3 MeV is very close to the 
frequency of oscillations obtained with the NL3 parameter set: 14.1 MeV. As in the case 
of ^°Zr, the calculated excitation energy for the NLl parameter set (-ft^nm = 211.7 MeV), 
is approximately 1 MeV lower than the average experimental value. 

The effective compression modulus (^4|) '^^ ^osp]^ jg fiigpiayed in Fig. 5 as a 
function of -R'nm- The observed behavior of is almost identical to the constrained 
compressibility calculated in the relativistic mean-field model using the generator coor- 
dinate method 0. The deviation from the almost linear behavior, which is observed for 
the HS linear parameter set, is slightly more pronounced in the present time-dependent 
calculation. 

The transition densities that correspond to the main peaks calculated for the six 
effective interactions are displayed in Figs. 6a and 6b. In the left columns we plot the 
dynamical transition densities, calculated as Fourier transforms of the time-dependent 
baryon densities. On the right hand sides the scaling densities are shown 

pf)(r) = 3pr^ + r^pf)(r), (25) 

where Po''"^ denotes the ground-state vector density for protons (neutrons). In the scaling 
model the transition densities follow from a simple radial scaling of the ground-state 
density, both the central density and the surface thickness vary. We notice that the 
scaling transition densities are almost identical, with a possible exception for LI, for 
all effective interactions. They do not provide any information about the dynamics of 
isoscalar monopole vibrations. On the other hand, the radial shape of the transition 
densities that result from time-dependent calculations depends very much on the value 
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of the nuclear matter compression modulus. As we have already seen for ^°Zr, starting 
with NL2 (-ft'nm = 399.2 MeV), a minimum develops in the center of the nucleus. It also 
appears that the oscillations of proton and neutron densities are out of phase at r = 0. 
Thus not only does the frequency of oscillations increase with -ftTmn, but also the radial 
dependence of the density oscillations changes dramatically. 

In order to understand better this behavior, we separate in Fig. 7 the volume and 
surface contributions to the transition densities. Following the procedure of Ref. [Q] we 
define a velocity field associated with the collective motion 

u(r) = -—^rr'^pj.{r\Eo)dr' (26) 
r-^Po{rj Jo 

where po is the ground-state density, and prif^', Eq) denotes the transition densities shown 
in Fig. 6. The transition density is separated into two components 

p??^(r,i?o) = PoVu = Po-,Ur^u) (27) 

pr\r,Eo) = uVpo = u^. (28) 

ar 

The resulting volume and surface transition densities are shown in Fig. 7 for all six 
effective interactions. We notice that the surface contribution does not depend much on 
the parameter set used, that is, on the nuclear matter incompressibility. The volume 
transition density, as one would expect, is very sensitive to the value of -ft'nm- A very 
interesting phenomenon is the formation of standing waves in the bulk. It starts already 
for NL2, but is clearly observed for HS and LI. 

The effective interactions NLl and NL3 seem to produce GMR excitation ener- 
gies which are quite close to the experimental values. We have therefore calculated, for 
these two parameter sets, the isoscalar giant monopole resonances in a number of doubly 
closed-shell nuclei: ^^Ni, ^°°'^^^Sn, ^^^Zr, ^"^^Gd. The results, together with those already 
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discussed, are shown in Fig. 8. The energies of giant monopole states are determined from 
the Fourier spectra of the time-dependent monopole moments, and are displayed as func- 
tion of the mass number. We notice that the NLl excitation energies are systematically 
lower, but that otherwise the two effective interactions produce very similar dependence 
on the mass number. The empirical curve ~ 80 A^^^^ MeV is also included in the fig- 
ure, and it follows very closely the excitation energies calculated with the NL3 parameter 
set. 

Experimental data on isovector giant monopole resonances are much less known. 
The systematics of excitation energies will not, in general, depend on the nuclear matter 
compression modulus. More likely, energies will depend on the coefficient of asymmetry 
energy. Isovector excitations are therefore outside the main topic of the present study. 
Nevertheless we have calculated the eigenfrequencies of isovector monopole modes, and 
compared them with available data on energies of giant monopole resonances. The time- 
dependent isovector monopole moments 

2 2 

< r > — < r > 

p n 

and the corresponding Fourier spectra for ^°^Pb are displayed in Fig. 9. Compared to the 
isoscalar vibrations (Fig. 4), we notice that the isovector oscillations are more damped and 
the anharmonicities are more pronounced. Consequently, the Fourier spectra are more 
fragmented. In each figure only the energy of the main peak is displayed. The four non- 
linear parameter sets produce similar Fourier spectra, with most of the strength between 
25 MeV and 30 MeV. The fragmentation of the Fourier spectra is more pronounced for 
the linear effective interactions HS and LI, with a considerable amount of strength shifted 
in the energy region 30 — 40 MeV. The results of time-dependent calculations should be 



compared with the experimental value for the IV GMR in ^°^Pb: 26 ± 3 MeV We 
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have also calculated the time-dependent isovector monopole moments for ^°Ca and ^°Zr. 
Although the oscillations are more anharmonic as compared to ^^'^Pb, the corresponding 
Fourier spectra are generally in agreement with the experimental data on the IV GMR: 
31.1±2.2 MeV for ^°Ca and 28.5±2.6 MeV for ^oZr [p. 

4 Constrained relativistic mean-field calculations 

A very useful method for description of excited states in nuclei is provided by con- 
strained mean-field calculations. In the framework of the relativistic mean-field theory, 
constrained calculations were performed for the ground state of ^^Mg, using two quadratic 
constraints p2|. By analyzing the curvature of the energy surface near the ground state, 



in Ref. constrained RMF calculations were used to determine the dependence of the 
excitation energy of giant monopole states on the nuclear compressibility. In Ref. this 
approach has been generalized by applying the generator coordinate method (GCM) to 
the RMF. The GCM takes into account correlations produced by collective motion of the 
nucleons. 

In the present study we further extend the method of Ref. [^, by using a more 
general ansatz for the generating functions of the GCM. Specifically, in Ref. [^j Slater 
determinants of nucleon single-particle spinors, resulting from constrained RMF calcula- 
tions, were used as generating functions. Here the bosonic sector is explicitly taken into 
account. The generating functions will be defined as direct products of Slater determi- 
nants built from single-particle spinors and the coherent states that represent the meson 
fields. Furthermore, the investigation is extended to include the isovector (T = 1) giant 
monopole states. 
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The GCM A-particle trial wave- function \l/GCM(ri5 • • • , ^a) is written in the form 

^GCM(ri, . . . , r^) = J f{q) ^(ri, . . . , r^; g) , (29) 

where \l/(ri, . . . , r^; g) are the "generating functions", and f{q) is the "generator" or 
"weight function" that depends on the "generator coordinate" q. From the variation of 
the energy of the system with respect to f{q), the //i//- VF/iee/er integral equation [^] for 
the weight function /(g) is derived 

J [n{q, q') - EM{q, q')] f{q') dq' = , (30) 

where 

n{q,q') = {^iq)\H\^{q')), (31) 

and 

Xiq,q') = {^iq)\^iq')), (32) 

are the Hamiltonian and normalization kernels, respectively. The solutions of the Hill- 
Wheeler equation are the discrete eigenvalues En and eigenfunctions fn{q), which deter- 
mine the nuclear wave-functions for the ground and excited states (for details see Ref. 



The relativistic extension of the GCM uses the self-consistent solutions of the equa- 
tions of motion of the Lagrangian density (|l]). The generating functions are defined as 
products of Slater determinants $(g), built from constrained RMF solutions for the Dirac 
single-nucleon spinors ipii^r; g), and the coherent states of the cr, u, and p meson field and 
the electromagnetic field 

l^'(g)) = |<l>(g)) ® |a(g)) ® \uj^iq)) ® \p^iq)) ® |A^(g)) . (33) 

The general expression for the boson coherent state reads pB| : 



.) = C 2 exp 
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|0) , (34) 



where 0^ denotes the a, uj^, and fields, ^/^(k) is the Fourier-transform of the field 

</)^(k) = uj{k) I ^^V0^(r)e*'^^ (35) 



u;{k) = ko = v/k2 + m2^ , (36) 

and 

(27r)32cu(A;) (27r)4^ ^ ^ "^"^ ^ ' ^ ' 

a"'"(k) is the creation operator of a boson with momentum k and C is the normalization 
constant. Onishi's formulas are used in the calculation of the integral kernels (|3T|) 
and ( p2D in the solution of the Hill- Wheeler equation. The normalization kernel can be 
written as product of fermionic and bosonic factors 

X{q,q') = XB{q,q'Wu{q,q'), (38) 

where 

^^J,{q,q') = mmq')) = det{iV,,(g, g')}, (39) 

with the matrix 

N^J{q,q') = iUqMiq')), (40) 

where I'lpiiq)) denotes the single-particle spinors. The boson factor 

^fM{q,q') = {cr{q)W{q')) Mq)Mq')) {pM)\pM')) {M<i)\M<i')) (41) 
can be calculated from the overlap kernel 

M)\<t^M')) = exp(-^|d3fc|0^(k;g)-0^(k;g')p). (42) 



The Hamiltonian kernel Ti.{q, q') reads 

n{q,q') = Af{q,q') J n{r- q, q') dr, (43) 
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where, for the case of spherical symmetry, the energy density kernel 7i(r; g, q') is given as 

A 

^RMF(r; g, q') = J2 ^jl' "Ptir; q){-iaV}ijj{r; q') + Mp,(r; g, q') 

+ \9^Ps{r; q, q')(y{r\ q, q') + ^gu^pv{r; q, q')uj^{r; q, q) 

1 1 

+ ^QpPzij] q, q')p°{r; q, q') + -epp(r; q, q')A^{r; q, q') 

+ ^((Va(r; g, g'))' + m>2(^; q')) + l^^^a' + ^(^s^^' 
-\{{Vuj\T-q,q')f + ml{u:\r-q,q')f) 
~\{{^p\r-q,q')f + ml{p\r-q,q')f) 

-i(VA°(r;g,g'))^ (44) 

with the density matrices 

A 

Ps{r;q,q') = J] iV^i ^,(r; g)^,(r; g') (45) 
P.(r;g,gO = ^ iy,:iV+(r; g)V^,(r; g') (46) 

A 

P3(r;g,g') = J] iV^i V;+(r; g)r3^,(r; g') (47) 
Pp(r;g,gO = iV^i ^+(r; g)il^^,.(r; g') • (48) 

i,j=l ^ 

The densities are sources for the Klein-Gordon equations. Their solutions are the meson 
fields that appear in Eq. (^) 

0/.('^;g,g') = ^(0/.(r;g) + 0;.(r;gO. (49) 

5 Isoscalar and isovector Giant Monopole Reso- 
nances 

Solutions of constrained RMF calculations have fixed expectation value Q of some ob- 
servable Q. The constraining operator Q is coupled to the Dirac Hamiltonian through a 
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Lagrange multiplier q 



{hD- qQ)ipi{Y]q) = eiiJi{r;q), 



(50) 



with the Dirac Hamiltonian 



hi) 



icxV + /?((M + g^a) + g^^Yu^ + gpi^fp^ + e j 




2 



(51) 



The constrained Dirac equation and the corresponding Klein-Gordon equations for the 
fields are solved self-consistently. In relativistic GCM calculations the Lagrange multiplier 
q is chosen as the generator coordinate. The constraint operator reads 

Q ^ = {2Urrl T = 0, 1 (52) 
1=1 

and it relates the generator coordinate to nuclear rms radii 



Using a Hartree code in coordinate space, constrained RMF calculations for the spherical 
symmetric case are performed for a number of values of the Lagrange multipliers q. The 
allowed values of q are limited by the requirement that the coupled system of Dirac and 
Klein-Gordon equations converges to a self-consistent solution. 

From a set of constrained RMF solutions (single-particle spinors of occupied states, 
meson and electromagnetic fields) , the normalization (^) and energy (|3T|) kernels are cal- 
culated. Solutions of the associated integral Hill- Wheeler equation determine the nuclear 
ground and excited states. The method of "symmetrical orthogonalization" [p6| , p7| is 
applied to the generalized eigenvalue problem. The Hamiltonian kernel is projected onto 
the regular (collective) subspace of the norm kernel. A series of collective eigenenergies of 
the constrained system and the corresponding nuclear wave-functions are obtained, which 




(53) 
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can be interpreted as the correlated ground-state solution and excited monopole states. 
Furthermore, nuclear densities in the ground-state and first excited state are calculated, 
as well as the transition densities 

Pfi{r) = J dqdq' ff{q)fi{q')Af{q,q')p,{r;q,q') , (54) 

where / and i denote the resonance and the ground-state, respectively. Using the four non- 
linear effective interactions NLl, NL3, NL-SH and NL2, calculations have been performed 
for a number of spherical nuclei. 

In Table 2 the isoscalar monopole excitation energies are displayed. For the NLl 
parameter set we also include in the second row the excitation energies obtained with 
constrained GCM calculations from Ref. |^. It appears that the inclusion of meson 
coherent states in the generating functions lowers the calculated excitation energies by 
approximately 0.5 MeV. Compared to the results of time-dependent calculations of Sec. 3, 
the energies in Table 2 are lower by more than 1 MeV, and is some cases the difference 
is larger than 2 MeV. The reason is that in constrained GCM calculations one implicitly 
assumes that the motion is adiabatic. In the time-dependent calculations on the other 
hand, no assumption about the nature of the mode of vibration is made. The comparison 
between constrained and time-dependent mean-field calculations strongly emphasizes the 
fact that we have to understand how the nucleus vibrates in order to relate the excitation 
energies of monopole states to the nuclear matter compression modulus. 

In Fig. 10 the calculated energies of giant monopole states are plotted as functions 
of A~^^^. An average linear dependence is observed for all four parameter sets. The devia- 
tions from a pure linear dependence are relatively small for NLl and NL3. They are more 
pronounced for the effective interactions NL-SH and NL2, for which the nuclear matter 
compression modulus K^ra > 300 MeV. The effective nuclear compression moduli Ka (|2^) 
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of ^^O, ^°Ca, ^°Zr and ^osp]^ g^j,g shown in Fig. 11 as functions of K^^. The GCM results 
for Eq and < >o are used in (24). We did not calculate from the curvature of the 
energy surface, and then Eq from (|2^) , as was done in Refs. [23, |5|. The relation (24) is 
used to obtain an estimate for the values that result from GCM calculation, in which 
the dynamics of the system is described by the non diagonal matrix elements of the energy 
kernel. Due to the correlations that these matrix elements produce, the calculated values 
of Ka are somewhat lower than the corresponding results of Ref. |Q. For ^osp]^ ^j-^g values 
of Ka can be compared with results of time-dependent calculations (see Fig. 5). The con- 
strained GCM transition densities (|MD for ^osp]-, ^re shown in Fig. 12. Similar in shape 
to those calculated in the scahng model (Fig. 6), the constrained transition densities dis- 
play a radial behavior which does not significantly depend on the effective interaction. In 
contrast to the results of time-dependent calculation, the dynamics of isoscalar monopole 
vibrations which is described by the constrained transition densities does not depend on 
the effective compression modulus. Again, this is due to the assumption of adiabatic 
motion, which is inherent in the constrained GCM approach. 

In Table 3 we display the excitation energies of isovector monopole states in '^^Ca, 
^°Zr and ^'^'^Pb, and compare the theoretical values with experimental data on isovector 



giant monopole resonances |21|. The calculations have been performed for the four non- 
linear parameter sets, and in the second column the corresponding asymmetry energy 
coefficients agym are included. We notice that the calculated excitation energies, like the 
results of time-dependent calculation, do not depend very much on the effective inter- 
action. In particular, no simple relation can be established between the coefficient of 
asymmetry energy and the excitation energies of isovector states. The calculated values 
for ^°Ca and ^''Zr are close to the experimental excitation energies. For ^''^Pb, however, 
constrained GCM calculations produce results which are almost 10 MeV lower that the 
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experimental IV GMR. This is in contrast with the results of time-dependent calculation 
which center around 28 MeV, in agreement with experimental data. This result might 
indicate that in a heavy nucleus such as ^osp]-, ^j^g assumption of adiabatic motion is not 
justified for the isovector mode. On the other hand, the constrained GCM results might 
also indicate that some strength should be expected around 17 MeV for the isovector 
monopole mode in ^'^^Pb. The corresponding time-dependent results (Fig. 9) display 
some strength below the main peak in the Fourier power spectra, but only very little 
below 20 MeV. Finally, the isovector transition densities for ^^^pio are shown Fig. 13. As 
we have already seen for the isoscalar mode, the transition densities depend very little on 
the effective interactions. 

6 Summary 

Relativistic mean field theory has been used to investigate the monopole eigenmodes of 
a number of spherical closed shell nuclei, from ^^O to the heavy nucleus ^°^Pb. Two 
microscopic methods have been applied: 

a) Time-dependent relativistic mean-field calculations. Starting from the self- 
consistent mean-field solution for the ground state and a set of appropriate initial 
conditions, the full set of time-dependent RMF equations is solved and dynamical 
variables are analyzed. Collective isoscalar and isovector monopole oscillations are 
studied, which correspond to giant resonances in spherical nuclei. Excitation en- 
ergies and the structure of eigenmodes are determined from a Fourier analysis of 
dynamical monopole moments and densities. 

b) The generator coordinate method, with generating functions that are solutions of 
constrained RMF calculations, is used to calculate excitation energies and transition 
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densities of giant monopole states. Constrained RMF calculations, with (T=0 
and T— 1) as the constraint operator, determine the stationary wave-functions with 
different values of nuclear rms radii. Using the generator coordinate method, the 
RMF-Hamiltonian is diagonalized in the basis of these solutions. The lowest eigen- 
mode corresponds to the correlated ground state. The solution for the first excited 
state describes the giant monopole resonance. 

Calculations have been performed for six different parameter sets of the Lagrangian. 
Two of them correspond to older effective forces, which do not include nonlinear self- 
interaction of the (T-field. The others are more modern parametrizations, which have been 
used with great success in recent years for a unified description of properties of nuclei over 
the periodic table (binding energies, radii, isotopic shifts, deformation parameters, mo- 
ments of inertia in superdeformed nuclei, nuclear halos at the neutron drip line). The six 
parameter sets differ essentially in the deduced value for the incompressibility of nuclear 
matter. Therefore, from the energy spectra and transition densities calculated with these 
effective forces, it has been possible to study the connection between the incompressibility 
of nuclear matter and the breathing mode energy of spherical nuclei. 

For the isoscalar mode we have found an almost linear relation between the excitation 
energy of the monopole resonance and the nuclear matter compression modulus. GCM 
calculations based on constrained solutions of the RMF equations yield, in general, shghtly 
lower excitations energies as compared with results of time-dependent calculations. This 
can be understood from the assumption of adiabatic motion which is implicitly contained 
in constrained RMF calculations. No assumption about the nature of the mode is made in 
the time-dependent model. With both methods, the effective interactions NLl and NL3 
produce GMR excitation energies which are close to the experimental data. Based on the 
very precise experimental value for the isoscalar GMR energy in the heavy nucleus ^°^Pb, 
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which is ratlicr well reproduced by the set NL3, we have derived a value of approximately 
270 MeV for the incompressibility of nuclear matter. For a series of doubly closed-shell 
nuclei, the effective interaction NL3 reproduces the empirical result ~ 80 A"^^^ MeV. 

Transition densities provide valuable information on the structure of the different 
modes. It turns out that there are essential differences between the two methods. The 
constrained GCM transition densities are similar to those of the simple scahng model, and 
display a radial behavior which does not significantly depend on the effective force. This 
is not the case for the time-dependent RMF calculations. Here the radial shape of the 
transition densities depends very much on the value of the nuclear matter compression 
modulus. Only for the parameter sets NLl, NL3, and NLSH, the transition densities 
display a radial shape characteristic for the breathing mode. For the other sets we find 
essential deviations in the nuclear interior, sometimes even standing waves. Only the 
surface component of the transition density is rather independent on the effective force. 

As expected, for the isovector mode we have found no relation between the ex- 
citation energy of the giant resonance and the calculated incompressibility of nuclear 
matter. More damping and more fragmentation in the spectra is observed, compared to 
the isoscalar case. Excitation energies that result from time-dependent calculations are 
in good agreement with available experimental data. The constrained RMF results agree 
with experiment only for light and medium heavy nuclei. For ^°^Pb we find results which 
are almost 10 MeV lower that the experimental T=l GMR. This might indicate that in 
heavy nuclei the assumption of adiabatic motion is not justified for the isovector mode. 

Summarizing our investigations in the framework of relativistic mean field theory, we 
conclude that the nuclear matter compression modulus Knm ~ 270 MeV is in reasonable 
agreement with the available data on spherical nuclei. 
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Figure Captions 

Fig. la Time-dependent isoscalar monopole moments < > {t) and tlie corresponding 
Fourier power spectra for ^°Ca. The parameter sets are NLl, NL3 and NL-SH. 

Fig. lb Same as in Fig. la, but for the effective interactions NL2, HS and LI. 

Fig. 2a Time-dependent isoscalar monopole moments < > (t) and the corresponding 
Fourier power spectra for ^°Zr. The parameter sets are NLl, NL3 and NL-SH. 

Fig. 2b Same as in Fig. 2a, but for the effective interactions NL2, HS and LI. 

Fig. 3 Transition densities for the isoscalar monopole states in ^°Zr. 

Fig. 4a Time-dependent isoscalar monopole moments < > (t) and the corresponding 
Fourier power spectra for ^^'^Pb. The parameter sets are NLl, NL3 and NL-SH. 

Fig. 4b Same as in Fig. 4a, but for the effective interactions NL2, HS and LI. 



Fig. 5 The compression modulus Ka ( plf ) of ^°^Pb as function of the nuclear matter 



compression modulus -ft'nm- The excitation energies result from time-dependent 
relativistic mean-field calculations. 

Fig. 6a Transition densities calculated as Fourier transforms of time-dependent baryon 
densities (left column), and scaling transition densities (right column), for the 
isoscalar monopole states in ^^'^Pb. The parameter sets are NLl, NL3 and NL- 
SH. Solid lines denote proton densities, and dashed lines correspond to neutron 
transition densities. 

Fig. 6b Same as in Fig. 6a, but for the effective interactions NL2, HS and LI. 
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Fig. 7 Volume and surface transition densities for the isoscalar monopole states in ^°®Pb. 
Solid lines correspond to proton densities, and dashed lines denote neutron transition 
densities. 

Fig. 8 Excitation energies of isoscalar giant monopole resonances in doubly closed-shell 
nuclei as function of the mass number. The effective interactions are: NLl (squares) 
and NL3 (circles). The solid curve corresponds to the empirical relation ^ 80 A~^^'^ 
MeV. 

Fig. 9a Time-dependent isovector monopole moments < > — < r"^ > and the corre- 
sponding Fourier power spectra for ^''^Pb. The parameter sets are NLl, NL3 and 
NL-SH. 

Fig. 9b Same as in Fig. 9a, but for the effective interactions NL2, HS and LI. 

Fig. 10 Constrained relativistic GCM excitation energies of isoscalar monopole states in 
doubly closed-shell nuclei as function of mass number. The effective interactions are 
NLl, NL3, NL-SH and NL2. 

Fig. 11 The effective nuclear compression modulus Ka ( p4D of ^^O, ^''Ca, ^°Zr and ^'^'^Pb 
as function of K^^. The excitation energies of monopole states are calculated with 
the constrained relativistic GCM. 

Fig. 12 Constrained GCM transition densities for the isoscalar monopole states in ^°^Pb. 

Fig. 13 Constrained GCM transition densities for the isovector monopole states in ^°^Pb. 
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Table 1: Parameter sets for the effective Lagrangian. 





NLl 


NL3 


NL-SH 


NL2 


HS 


LI 


m [MeV] 


938.0 


939.0 


939.0 


938.0 


939.0 


938.0 


[MeV] 


492.25 


508.194 


526.0592 


504.89 


520.0 


550.0 


[MeV] 


795.335 


782.501 


783.0 


780.0 


783.0 


783.0 


THp [MeV] 


763.0 


763.0 


763.0 


763.0 


770.0 




9a 


10.138 


10.217 


10.44355 


9.111 


10.47 


10.30 


9uj 


13.285 


12.868 


12.9451 


11.493 


13.80 


12.60 


9p 


4.975 


4.474 


4.3828 


5.507 


4.04 




92 [im.-^] 


-12.172 


-10.431 


-6.9099 


-2.304 






93 


-36.265 


-28.885 


-15.8337 


13.783 






Knm [MeV] 


211.7 


271.8 


355.0 


399.2 


545.0 


626.3 



Table 2: Constrained GCM energies of isoscalar monopole states. The values of Knm and 
the excitation energies are in MeV. 







K 


160 


40Ca 


48Ca 


90Zr 


208pb 


1 


NL-1 


211.7 


20.2 


16.6 


15.9 


14.1 


11.0 


2 


NL-1 i 


211.7 


20.6 


17.1 




14.7 


11.7 


3 


NL-3 


271.8 


22.6 


19.6 


18.9 


16.9 


13.0 


4 


NL-SH 


355.0 


25.0 


22.0 


21.5 


19.5 


15.0 


5 


NL-2 


399.2 


27.1 


24.4 


23.0 


21.9 


16.6 



Table 3: Constrained GCM energies of isovector monopole states. The values of agym and 
the excitation energies are in MeV. 









40Ca 


902r 


208pj^ 


NL-1 




43.5 


29.0 


26..3 


16.5 


NL-3 




37.4 


28.6 


27.4 


18.0 


NL-SH 


36.1 


28.5 


27.9 


18.4 


NL-2 




43.9 


30.3 


28.8 


16.9 


exp. 


21| 




31.1±2.2 


28.5±2.6 


26.0±3.0 
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